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' Abstract 

o : 

■ The field algebra of the minimal models of W-algebras is amenable to a 

very simple description as a polynomial algebra generated by few elementary 

£S) ' fields, corresponding to order parameters. Using this description, the complete 

Landau-Ginzburg lagrangians for these models are obtained. Perturbing these la- 

43 grangians we can explore their phase diagrams, which correspond to multicritical 

points with D n symmetry. In particular, it is shown that there is a perturbation 
for which the phase structure coincides with that of the IRF models of Jimbo et 

43 ' al. 

> : 

•i-H . 

1 Introduction 

The classification of two-dimensional conformal field theories (2dCFT) is a flourishing 
branch of mathematical physics. It can be applied in essentially two different directions. 
One is string theory; the other is phase transitions in two dimensions and, hopefully, 
also higher dimensions. With regards to this second application, the information pro- 
vided by 2dCFT may seem excessive. In phase transition physics, before addressing 
the question of critical exponents and, therefore, dimensions of fields, one is usually 
more interested in determining the phase diagram. For this purpose, the essential 
information consists of the order parameters (with their symmetry) and the relevant 
fields formed out of them. We are used to seeing this information coming from the 
Landau potential, or its field-theory version, the Landau-Ginzburg lagrangian. This 
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lagrangian is not a datum of 2dCFT but it can be obtained from it, as was shown by 
A.B. Zamolodchikov |l| in the simple case of the unitary series of minimal models. All 
these models have only one order parameter and Z2 symmetry. In order to describe 
phase transitions with higher symmetry, one is compelled to consider 2dCFT with an 
extended algebra. Among them, the minimal models of the W-algebras are the 
natural generalization of the Virasoro minimal models. 

The Zamolodchikov's procedure to associate lagrangians to 2dCFT begins by giving 
a polynomial structure to the algebra of primary fields. He found that all the relevant 
primary fields of a minimal model can be expressed as composite (powers) of the most 
relevant one (elementary field), which plays the role of order parameter. This is a 
necessary condition for the existence of a lagrangian and it is actually sufficient to 
guess it. Nevertheless, there is a direct way to obtain it, relying on the fact that the 
composite field next to the most relevant field has to be identified with the descendant 
of the elementary field, giving therefore the equation of motion. In this paper, we will 
generalize this procedure to the minimal models of the W-algebras. An attempt in this 
direction was already made by I. Koh and S. Yang [HI, but they did not obtain the 
complete structure of the algebra of primary fields and their lagrangians missed some 
important terms. 

The Theory of Catastrophes is relevant in the context of Landau potentials, as a 
mathematical framework in which many intuitive concepts of phase transitions can be 
precisely formulated. Its classificatory power has already found wide application in 2d 
N = 2 superconformal theories. However, it is still a rather specialized tool and we 
should not assume that the reader is familiar with it. Therefore, we will put it aside, 
though making some remarks in footnotes. 

The paper is divided in three parts. We begin with the simplest case, namely 
the W(3) algebra minimal models. In the first part, we explicitly identify all the rel- 
evant fields as powers of two elementary ones and therefore endow the algebra with 
a polynomial structure. The equations of motion and the lagrangian ensue from the 
identification of further powers of the elementary fields. In the second part, the phase 
structure that results from these potentials is briefly analyzed. The third part is de- 
voted to the general W-algebra models. To deal with them one has to take into account 
some new features that already appear for . The first model of the W^) series has 
central charge c = 1 and belongs to the Ashkin- Teller critical line |4| . For its impor- 
tance we take it as an example to study those new features. We finally discuss the 
phase structure of the general models and their relation with the solvable statistical 
models defined by the japanese group || (JMO models). 
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2 Structure of the field algebra of the W{%\ models 

The minimal models of the W(s) algebra are constructed in analogy with the current 
Virasoro minimal models, but using instead a two-component free field ||. The proce- 
dure relies on the Dotsenko-Fateev screening charge method. There are four screening 
charges; the primary fields are therefore labelled by a 2 x 2 matrix 




Every unitary model is still determined by one integer p and noted W%\ ■ This model 
contains p(p — l) 2 (p — 2)/12 independent primary fields and has D 3 symmetry. The 
fusion rules are derived from the neutrality condition including definite numbers of 
screening charges, 

ck( 3 ) = a(i) + a (2 ) + iVa+ei + N'a + e 2 + Ma_ei + M'a_e 2 , (1) 

with 

m(i)]a_)«i + 

m / (i) ]a_)w 2 . (2) 

Here e, and are the positive roots and fundamental weights of 577(3) (A 2 algebra), 
respectively. Note that Eq. ([!]) expresses the Clebsch-Gordan decomposition of the 
tensor product of two S77(3) ® SU(3) representations with highest weights — c*(i) and 
— a (2). The solution for , n', 3 \ , , m'/ 3 \ is 

n(3) = 7l(i) + 7l (2 ) - 1 - 2JV + JV' 

n{ 3) = n ' {1) +n{ 2) -l-2N' + N (3) 
and similar equations for m(3),m', 3 y 

2.1 Elementary and composite fields 

The most relevant fields are 

a = $ (? J)' * = 9 Q 2)' 

which represent the /J3 spin density and its conjugate |Q. They are the obvious can- 
didates for elementary fields. It is straightforward to prove that a certain number of 
the next most relevant fields are obtained from them by the Zamolodchikov's method 



7i(i) 



([1 - n (i) ]a + + [1 
([l-n' (i) ]a+ + [l 
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with the fusion rule ([!]). Besides, they arrange themselves in a triangular structure 
corresponding to the SU(3) lattice of dominant weights. Before proceeding further, let 
us recall Zamolodchikov's field identifications in more detail. The primary field content 
of the Virasoro minimal model M p can be visualized on a grid of dimension px (p—1) 
symmetric with respect to the center. The elementary field is placed at (2, 2) and its 
powers are placed on the main diagonal, up to (p — l,p — 1). The next power is at 
(p,p — 2) or equivalently at (1, 2). This is the lowest end of the second diagonal, which 
contains the remaining powers, up to the 2p — 4th, corresponding to the least relevant 
field. 

In the present case, the conformal grid is four-dimensional and the equivalent of 
the "main diagonal" is a two-dimensional triangular section of it. To be precise, this 
"main diagonal" exactly corresponds to the SU(3) lattice of dominant weights of level 
p — 3, as it is proved by the identifications 

with 0<k + l<p-3, obtained from (g) for the trivial N = N' = M = M' = case. 
The next power, k + I = p — 2, leads to essentially three different choices. For 1 = 0, 
k = p — 2, the appropriate choice is N = 1, N' = M = M' = 0, and 

This choice holds for the subsequent powers filling the triangle subtended from it up 
to the 2p — 6th power, giving 

W.^J HI), (6, 

with 0<k + l<p — 4. There are similar identifications for a p ~ 2 and the triangle 
conjugate to the one mentioned above, with N = 0, N' = 1, M = M' = 0. We are left 
with the middle triangle, which requires N = 1, N' = 1, M = M' = 0. A convenient 
form of the ensuing identifications is 

a p-k-2- k+ i-i = J k+iy (7) 

with 1 < k, I and k + I < p — 2. In particular, for k = 1, 1 = p — 3, 
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we get the least relevant field of the thermal subalgebra (containing the fields with 
n = n' and m = m'). This subalgebra is actually generated by its most relevant field, 
the energy density 

'2 2 




6 — (7(7 

Its succesive powers span the subalgebra up to e p ~ 3 , identified above. This field will 
be important in the sequel. Its conformal dimension, as obtained with the general 
formulas ||, is (p — 2)/(p + 1). 

The three triangles containing the powers of the elementary field ranging from the 
p— 2nd to the 2p— 6th form a two-dimensional section of the conformal grid equivalent to 
the second diagonal of the Virasoro case. The field identifications above are consistent 
across the borders, although in a nontrivial way. The fields on the lower side of the 
upper triangle are 

= iili ili)- < 9 > 

They must me compared to the products of sigma with the fields on the upper side of 
the middle triangle 

' 1 2 N 



There is exact coincidince if the fusion with a is performed according to equations 
with .\/ 1..V .V M' 0. 




We know that in the Virasoro case, the main and second diagonals of the conformal 
grid contain all the relevant fields of the model. We have here a similar property ; 
namely, the big triangle with all the powers up to the 2p — 6th contains all the relevant 
fields of the W% model, as can be checked with the dimension formula. 



2.2 Field equations and Landau-Ginzburg lagrangians 

The fields corresponding to the 2p — 6th power are specially important. The most 
relevant is e p ~ 3 . Pursuing the analogy with the Virasoro case, we expect that further 
products with a or a must include the lowest dimension descendants, namely da or 
da, and therefore originate just field equations. However, in contrast with the Virasoro 
case, other fields already identified appear before these descendants in the regular part 
of the operator product expansion (OPE); they have to be kept in the field equations. 
This new feature was already remarked for the p = 4 case in J7[. The OPE for e p_3 in 
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particular, can be written as 



1 I J \1 2 J Z A #- 3 
1 1 



+ A CT+ A,„_,-A, 4» + ,A,„_ 3 -1 d<T+---, (12) 



where we have shown the most singular term and the initial regular terms; numerical 
coefficients have been omitted. The field 0, yet to be determined, must have been 
already identified as a power of a and a, and needs to be such that the corresponding 
exponent of z is positive and smaller than 1 — A eP -3, that is to say 

4 

A. < 1 + A„ = ' 



3p(p + 1) 

We notice that the fusion rule Eq. ||| with M — 1, N — N' — M' = 0, gives 

ae^ 3 = $^ 3) = 2 3 P 1 3 ) = C7 ^ V ~ 2 ( 13 ) 
with conformal dimension 

J2 A = 3p(p-l) + 4 

ll 3j 3p(p + l) 1 J 

This field therefore satisfies the condition and can be taken as <j). 

With <f) identified as above, the OPE fll2|) yields the following field equation (Hence- 
forth, the symbol ~ will stand for equality up to numerical coefficients) 

d 2 <r ~a p - 4 a p - 2 + a(aa) p - 3 . (15) 

This equation leads to the lagrangian 

C ~ dada + a^&P- 1 + (aa) p - 2 + c.c. 

= dada + (aa) p - 4 (a 3 + a 3 ) + (aa) p - 2 . (16) 

The new composites appearing in it can be identified as 

{™y- 2 = 4l \), (17) 



a p-4-p-l = $ ( 



■1 2 

4 2 ' 



which are of course irrelevant fields. The first one is the least irrelevant thermal field 
e p ~ 2 of conformal dimension A eP -2 = 1 + (3/p). The other is slightly more irrelevant, 
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A = 1 + [3(p + 2) / p(p + 1)], and is essential to endow the lagrangian with the necessary 
D 3 symmetry. Fields with higher dimension must also appear in the lagrangian, for 
there are more field equations coming from the products of a (a) and other fields with 
the 2p — 6th power. Those additional fields have to be neutral, and correspond to 
monomial terms with 2p — 4th and 2p — 5th powers, 

{aay- 5 {a 6 + a 6 ),{aa) p - 7 {a 9 + a 9 ), 

etc. These terms are important for the potential to be well defined^ but, since they 
represent more irrelevant fields, a renormalization-group argument would tell us that 
their effect becomes negligible close to the multicritical point. This circumstance can 
be mimicked by assigning them small coefficients. Therefore, they are not expected 
to play any relevant role in the phase structure of the models and will be ignored 
henceforth. Nevertheless, the term (<to") p-4 (<7 3 + (x 3 ) in (|1^) needs to be kept, not only 
to enforce the D 3 symmetry but also to prevent the potential from being degenerate. 
The degeneracy of the potential consisting of just the (aa) p ~ 2 term (the one found in 
||) manifests itself as the possibility of obtaining under some perturbations an infinite 
number of minima (one or several circumferences). Summarizing, the lagrangian ([TB]) 
is the correct starting point for the study of the phase diagram in the next section. 

We would like to point out a peculiarity of the lowest-p models: Counting their 
primary fields one realizes that there are not enough to support the identifications 
above. For instance, Wh\ and WA\ have 6 and 20 primary fields, respectively, which 
cannot account for all the fields mentioned above. The six fields of the former are all 
relevant and identifiable with the composite fields up to second power of the elementary 
ones. However, we know that this model coincides with the nondiagonal modular 
invariant (A4, D4) of the Virasoro series [|, which has additional primary fields for 
higher powers of the elementary fields, namely, 

$(12,13) = <7 3 , $(13,12) = O 3 (19) 

$(i3) = e 2 = {oaf. (20) 

Therefore, these fields have to be W-secondaries. It is easy to find that indeed 

a 3 = W^e, a 3 = W^e (21) 
e 2 = W_i W_i e. (22) 

In conclusion, some of the composite fields of the low-p models turn out to be W- 
secondaries with the necessary symmetry properties. 

1 This fact is specially clear in the language of Catastrophe Theory, where the absence of those 
terms produces an indeterminate potential H. 
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3 Phase structure of the W($\ models 



The mean-field phase diagram is obtained analyzing the various configurations of min- 
ima produced upon perturbation of the potential part of the lagrangian by relevant 
fields. The topology of the mean-field phase diagram is identical to the actual phase 
diagram (this is the reason why the Landau potential suffices to study the phase struc- 
ture). A complete representation of the phase diagrams is however extremely compli- 
cated, as can be imagined from the difficulties already encountered for the simplest 
case, the three-state Potts model Wfcs, worked out in 0. Nevertheless, one is mostly 
interested in certain properties of the phase diagram, namely, the number of possible 
phases (stable states) and their topological interrelation, which only demand a limited 
knowledge of it. These properties are conveniently expressed by the state diagram, 
which we define as a set of points corresponding to the possible states (either stable 
or unstable) in a definite region of the phase diagram, and being linked whenever they 
can merge under a suitable perturbation^. This diagram provides a direct connection 
with statistical mechanics: A perturbed conformal model can be interpreted as the 
renormalizat ion-group universality class of lattice models with a state variable living 
on that diagram. Now we recall that the W? n \ are known to describe the critical be- 
havior of the JMO lattice models || (see also [fTOfl). These are interaction round a face 
(IRF) models with state variables defined on the dominant-weight lattice of some Lie 
algebra, as a generalization of the restricted solid on solid (RSOS) models of Andrews, 
Baxter and Forrester, the critical behavior of which is described by the Virasoro min- 
imal models. A motivation for the construction of these IRF models actually arose 
from the representation of the minimal W models as cosets (A n © A n , A n ) with level 
(/ — 1, 1), in which the A n dominant weights of level / play a natural role (recall the 
role of (U) as the Clebsch-Gordan decomposition rule). 

The relation between the state diagram of the model and the diagram of 
dominant fundamental weights of A 2 is apparent in (7|. It is to be expected that the 
state diagram of the W?$\ model coincides with the A 2 dominant-weight diagram of level 
p — 3. Therefore, we would be interested in obtaining that state diagram by perturbing 
the potential in a suitable way. Finding this perturbation in the general case turns 
out to be quite a nontrivial problem: Contrary to the Virasoro case, with one order 
parameter, for which is rather easy to discern the effect of each individual perturbation, 
in the case with several order parameters it is almost impossible to predict how the 
potential will deform under a particular perturbation. However, there are some rules, 
as we shall see. First of all, let us introduce real variables, by a = x + iy (a — x — iy), 

2 Hcnce, for phase diagrams derived from a potential the state diagram coincides with the Dynkin 
diagram of a level curve, as defined in || and used in this context in J?]]. However, we prefer here the 
name state diagram since we will be using Lie-algebra terms, among which Dynkin diagram already 
has a very standard meaning, not identical to that in ||. In fact, the state diagrams for the W( n ) 
models are related with weight diagrams of the underlying Lie algebra. 
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more suitable for representation of the potential, and simplified notation for the two 
independent D 3 invariants 

I = aa = x 2 + y 2 , 
h = cr 3 + a 3 = x 3 - 3xy 2 . 
The potential takes the form of a polynomial in Jo and 7i, 

V(x, y) = II 2 + 11% + P(I , 70, (23) 

where the symmetric perturbation V(Iq, I±) has degree 2p — 6 as a polynomial in x, y. 
If the potential contained only 7 , it would have full 0(2) symmetry and depend on 
just the radial variable r = y/x 2 + y 2 . Although this possibility is not admissible in the 
present context, it is convenient to compare with the case with one order parameter. 
Here the maximum unfolding deformation is achieved for non-null values of all the 
coupling constants in the symmetric perturbation. The modification introduced by 
Ii amounts to reduce the symmetry of this rotational invariant potential to D 3 . The 
structure is still the same along the three directions in the (x, y) plane for which Ii = 0, 
and must be qualitative similar all over. Therefore, one has to tune carefully all the 
coupling constants to produce minima configurations corresponding to the desired state 
diagram. This problem, topological in essence, cannot be formulated in any simple 
algebraic form. The simplest models can be succesfully handled by computer programs 
capable of representing algebraic surfaces, such as Mathematica. In this way, one tests 
different values of coupling constants to observe how the potential deforms. With 
some insight and a good deal of luck, one may arrive to a configuration recognizable as 
described by the appropriate A 2 dominant-weight diagram. Needless to say, the number 
of coupling constants increases rapidly with p and beyond p = 6 this trial-and-error 
method ceases to be feasible. 

Fortunately, there are more powerful topological arguments that can be directly 
applied to solve the inverse problem; namely, finding from the state diagram a potential 
with a minima configuration that fits it. Afterwards, this potential can be compared 
with the original one. In order to expose the method, let us introduce the concept of 
level curves for a potential; in our case, V(x, y) = c defines an algebraic curve at levelQ 
c. This curve is singular if there are extremal points of the potential (minima, maxima 
or saddle points) at that level. Now we consider the A 2 dominant-weight diagram of 
level I; we mark the middle points on each link and construct a singular curve joining 
these points, as in Fig. 1 Q. We would like to identify this topological curve with a 
level curve of the potential (^3|) for I = p — 3. The first thing we observe is that the 
curve is the product of a definite number of components (two in Fig. 1). We further 

3 Throughout the paper, we use the term level with two different meanings, namely, as the value of 
z = V(x, y) in the (x, y, z) coordinate representation providing the graphical notion of level curves or 
as the integer associated to a representation of a Lie algebra. 

4 This process is the inverse of obtaining the Dynkin diagram of a curve according to H. 
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notice that there are precisely (I + l)/2 trefoil shaped curves when I is odd and 1/2 
trefoil curves plus one of circular type when I is even. The level curve of the potential 
for the Wt\ model, 

V(x,y) = I 2 +h+wI = } (24) 

is the simplest algebraic curve with the trefoil shapef]; (observe that there are two 
different topological types depending on the sign of w). The simplest circular curve is, 
of course, the circle Iq — c = 0. Therefore, the simplest algebraic curve which represents 
the whole topological curve is given by the product 

(2+l)/2 

IJ (Io +h+ wJo) = I l +1 + I^Ii + • • • = 0, (25) 
i=i 

for /odd, and by 

1/2 

n(/o +h + WiI Q ){I -c)= P +1 + It^i + • • • = 0, (26) 

i=l 

for I even. We see that both reproduce the potential (|2"B| ) when I = p — 3. 

We should remark that the construction above yields values of the coupling con- 
stants (in fact, a range) for which the state diagram is the desired A2 dominant-weight 
diagram. However, there are certainly other values that produce the same diagram: All 
the level curves obtained by this method are singular and correspond to very special 
potentials with the saddle points at the same level; we can easily imagine deformations 
of the potential that move some saddle points off that level without altering its topol- 
ogy. There can actually be a simpler set of values, with some of them null. Moreover, 
those potentials are special in yet another sense: We have to contemplate the pres- 
ence of terms with higher powers of I\ in the non-perturbed potential, as we noted at 
the end of the first section. These terms, though, have small coefficients and will not 
change the topological type of the potential. In any case, that construction provides a 
privileged starting point to describe the phase diagram and can be used to probe the 
effect of other terms. We have done so with Mathematica, drawing various contour 
plots for the p = 5 and p = 6 models. A characteristic plot clearly exhibiting the state 
diagram appears in Fig. 2. 

4 The general W( n ) case 

The generalization of the Dotsenko-Fateev construction to W< n ) demands a free massless 
scalar field with n — 1 components. Therefore, the number of screening charges is 

5 This potential has been described in a previous paper 0. That (|24|) represents a trefoil can be 
understood without calculation when w = 0: I\ — gives just three straight lines intersecting at the 
origin and producing six sectors where V(x,y) takes negative and positive values alternatively. The 
effect of adding Ig is raising the potential for r large enough, hence closing the three negative sectors. 
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2(n— 1) and the primary fields are labelled by an equal number of integers. The unitary 
models WL\ contain [p!(p — l)!]/[n!(n — l)!(j» — n)\{j> — n — 1)!] spinless primary fields 
and they have D n symmetry |2|]. The neutrality condition to derive the fusion rules is 
the straightforward generalization of that in ([I]) 

n— 1 n— 1 

a(3) = a(l) + a(2) + a+ £ Afe + a_ £ M^, (27) 

i=l i=l 

with 

71-1 

a{j) = a({m(j)} | {rntO')}) = EQ 1 -^0')K + [1 - m^a-M- (28) 

i=l 

Here and uji are the positive roots and fundamental weights of SU{n) (algebra 
A n -x, respectively The fusion rule (1371) represents the Clebsch-Gordan decomposition 
for SU(n) <8> SU(n). To obtain the numerical fusion rules one is to substitute for the 
roots in terms of the weights, 

ei = 2uj 1 - u> 2 , 

e* = 2u>j - - u> i+ i, i = 2, ...,n-2, 
e n _! = 2o; n _! - w„_ 2 , (29) 

and solve for rij(3) and raj (3), 

m(3) = m(l)+n 1 (2)-l-2JV 1 + iV 2 , 

7ii(3) = n i (l)+n i (2)-l-2^ + iV i _i + iV i+1 , i = 2,...,n-2, 
n n _i(3) = n„_ 1 (l)+n„_ 1 (2)-l-2iV n _ 1 + iV n _ 2 (30) 

(There are similar equations for m*). 

The elementary fields are the n — 1 spin fields corresponding to the fundamental 
representations of SU(n), 

a k = $(1, . . . , 1,^, 1, . . . , 1 | 1, . . . , 1,^, 1, . . . , 1). 

k k 

They support the D n representations 

Ok — ► e 7 o-fe 

namely, the most relevant fields, o\ and o\ = a n -i, support the standard representation 
and the others the remaining two-dimensional representations; (if n is even, the central 
field forms a one- dimensional representation). The crucial fact is that all the other 
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relevant fields can be formed as powers of those n — 1 elementary ones and arranged on 
the SU(n) lattice of dominant weights of level 2(p — n). Showing this in the analogous 
way to the W(3) case, demands thinking of a 2(n — l)-dimensional conformal grid and 
various (n — l)-dimensional sections. We will restrict ourselves to n = 4 for simplicity, 
given that no essentially new properties arise for higher n. In this case we have three 
elementary fields, 

a 1 = $(1,0,0 | 1,0,0), 

2 = $(0,1,0 | 0,1,0), 

ox = $(0,0,1 | 0,0,1), 

the second one being real. The lattice of dominant weights is an isosceles pyramid 
(Fig. 3). The "main diagonal" contains fields up to level p — 4, which are identified in 
a straightforward way, from (|3~0"1) when all JVj = Mj = (i = 1,2,3), with the products 



of the elementary fields up to maximum total power p — 4, 

o\o\o™ = $(& + 1,1 + l,m + 1 | k + 1,1 + l,m + 1), k + l + m<p-A. (31) 

The identification of the next power, k + I + m = p — 3, demands that some of the iVj 
be non-null. For o\ and o\ we have choices analogous to those in the W@) case: 

• N% = 1 and the others null for 

a?" 3 = $(p - 4, 2, 1 | p - 2, 1, 1) = $(2, 1, 1 | 1, 1, 1). (32) 

• The conjugate equation when jV 3 = 1 and the others null. 

• The analogous middle triangle is also obtained by N\ = N% = 1 (the remaining 
null). 

However, the new component 02 introduces further possibilities, the simplest of which 
are 

of-% = 3,1,1 I p- 3,1,1) $(1,2,1 I 1,2,1) 

= $(p-4,l,l|p-3,2,l) = $(l,l,2|2,l,l), (33) 

with Ni = N2 = 1 (remaining null), and its conjugate. Another interesting one is 

cxf- 3 = $(l,p-3,l I l,p-3,l) $(1,2,1 I 1,2,1) 

= $(l,p-4,l|l,p-2,l) = $(l,2,l|l,l,l), (34) 

with Nx = N3 = 1, N2 = 2. All these choices hold for the respective domains of 
pyramidal shape that are subtended from the position of these fields in the big pyramid 
up to level 2p — 8 containing all the relevant fields. 
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An exhaustive description of the field identification domains and their matching 
conditions is cumbersome but not necessary to find the Landau-Ginzburg lagrangian: 
Like in the W(3) case, within the fields of highest level only the most relevant is needed. 
This field also belongs to the thermal subalgebra, now generated by two fields, 

ex = $(1,0,1 | 1,0,1) = (71^ (35) 

and 

e 2 = $(0,2,0|0,2,0)~(T 1 cTi + a 2 2 (36) 

(Recall that the symbol ~ means equal up to numerical coefficients). The seeked field 
is 

er 4 = (^ir 4 = $(M,l|2,l,2). (37) 
Upon product with the elementary field a±, we obtain the field equation 

c>Vi ~ ^r^rVs + Ma^y- 4 (38) 

(Note the similarity with (fRJ), although the monomial for the field <fr is now different, 
consistently with the present symmetry). Hence the lagrangian 

C ~ dcTtdat + (a^y-^al + a\)a 2 + {a^f^. (39) 

These are the crucial terms that determine the phase structure. However, like in the 
W(s) case, other terms corresponding to more irrelevant fields are also required for the 
potential to be well defined. In addition, we will consider again a symmetric pertur- 
bation, formed by lower degree terms corresponding to relevant fields. To simplify the 
aspect of all these terms, it is convenient again to regard the potential as a polynomial 
in the basic invariants. The invariants for the standard representation of D± and their 
expression in terms of real variables (a\ = x + iy) are 

J = o x o\ = x 2 + y 2 , 

I 1 = at + a A 1 =2(x i + y A -Qx 2 y 2 ). 

Now the presence of 02, forming a one-dimensional D4 representation^], produces new 
invariants 

I 2 = o\ = z 2 (z = cr 2 ) 

J 12 = (al + al)a 2 = 2{x 2 - y 2 )z 

That is to say, four basic invariants altogether. However, they are not independent, for 
they satisfy the relation 

4 = (h + 2/ )/ 2 . (40) 

6 The complete representation of D4 for the W($\ models is the direct sum Ai © E, in the notation 
of Landau-Lifshitz's textbooks. 
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This fact introduces some complications^ but it will suffice here to say that it amounts 
to the possibility of substituting for any power of lyi (higher than the first), which in 
consequence will not appear in the potential. Therefore, the potential takes the form 

V(x, y, z) = II 3 + F - 5 I 12 + V(I , h, h, I 12 ), (41) 

where the symmetric perturbation 

V(I , h, h, I12) = V (I , h, h) + InPx{h, h, h) 

has degree 2p — 8 as a polynomial in x, y, z. The irrelevant additional terms are 
Iq~ 5 Ii, Jg~ 4 / 2 , Io~ 5 Ii 5 Io~ 7 I^h, Io~ 6 Iih, etc., and will be ignored henceforth. 

4.1 Phase structure and its relation with IRF models 

In the study of phase structure we would like to generalize the method applied to 
the W(3) models. If it applies to the ones as well, we should expect to have 
configurations of minima that reproduce the succesive sets of level-/ dominant integral 
weights of A 3 . The simplest model, Whs, is again very illustrative. We will take 
advantage of the fact that there is an alternative description of this model: It is a 
special point on the line of 2dCFT defined on the orbifolds of a circle with variable 
radius. These 2dCFT, with c = 1 and D4 symmetry, correspond to critical Ashkin- 
Teller models with variable 4-spin coupling (0, [12|]). They have two twist fields with 



conformal dimension (1/16, 1/16), realizing the standard D4 representation (our x, y), 
and one field with dimension depending on the coupling constant K as 1/8ttK [Q, 
which realizes the non-trivial one-dimensional representation (in our model this field 
is z, with A 2 = 1/12). 

Specifying in ( |4"T| ) and after performing a tt/4 rotation on the xy- plane to write the 
cubic invariant in a more suitable form, the perturbed potential for the model can 
be written as 



V{x, y, z) = J 2 + I X2 + wl + w'I 2 = {x 2 + y 2 ) 2 + xyz + w{x 2 + y 2 ) + w'z 2 , (42) 



This potential is analogous to (|24j) except for the presence of two independent energy 
perturbations; the cubic term determines likewise its symmetry. The equation I 12 = 
represents the three coordinate planes dividing the space in eight sectors for which the 
potential takes minus and plus signs alternatively. Now the effect of adding Jq does not 
suffice to close a surface and produce minima, for there is no term z 4 , and the potential 
remains unbounded in this direction^. Nevertheless, the quadratic perturbation z 2 
bounds the potential, giving minima in pyramidal configuration. 



7 For a short account on Invariant Theory bearing on this problem, see |fll| . 

8 Note that the algebraic surface (x 2 + y + z 2 ) 2 + xyz = has a four-lobed tetrahedral shape and 
therefore posseses higher symmetry, namely, the full group of the tetrahedron Td (Landau-Lifshitz 
notation), isomorphic to the symmetric group S4. 
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Let us attempt to calculate the position of the minima of the full potential 
from the equations 

dV 

— = Axix 2 + y 2 ) +yz + 2wx = 0, (43) 
ox 

dV 

— = Ay(x 2 + y 2 ) +xz + 2wy = 0, (44) 
dy 

dV 

— = xy + 2w'z = 0. (45) 
oz 

The last equation is particularly simple to solve, giving 

z = -xy/2w', (46) 

which substituted back into the others yields the same equations as those coming from 
the two-variable potential with _D 4 symmetry 

V{x,y) = {x 2 + y 2 ) 2 + {2w')~ 1 x 2 y 2 + w{x 2 + y 2 ) 

= x * + y* + (2-^)x 2 y 2 + w(x 2 + y 2 ). (47) 
2w' 

We know that its extrema are placed on the vertices of a square. For the full 3-variable 
potential, the respective z obtained from ^ are identical in absolute value but negative 
or positive according to whether x and y have the same sign or not. Therefore, they 
form the vertices of an isosceles pyramid. However, it is not possible to check if the 
dimensions of this pyramid coincide with those of the A 3 dominant-weight lattice of 
level one (as drawn in Fig. 3). In fact, this question is nearly meaningless, given that 
there is no relation between z units and x or y units. Since the potential is at most 
quadratic in z, implying that this variable cannot affect critical behavior]], a healthy 
point of view is to disregard it altogether and consider just ( f47l) as the whole potential. 
This potential has already been associated to the critical line of the Ashkin- Teller 
model [14]. This line arises as a result of the existence of a marginal field, identified 
with the term x 2 y 2 , in that modelQ. Its coefficient a = 2 — (l/2w') is directly related 
to the coupling constant nK = (a/2) + 1, that is to say to the dimension of the field z. 
For a = we have two decoupled Ising models. When a = 2 the symmetry augments 
to 0(2), representing the Kosterlitz-Thouless point of the XY model. In principle, the 
parameter a can also take negative values as long as a > —2. For instance, the value 



9 This property is essential in Catastrophe theory, where a quadratic potential is the simplest Morse 
function, which can be added to any catastrophe to give an equivalent one. The particular equivalence 
to which we are led in our case discarding z is T244 ~ Xq of Arnold's classification jl3| . 

10 In Catastrophe Theory the potential (|47]) is called double-cusp catastrophe {Xg) and has aroused 
interest as the simplest catastrophe with a modal deformation flUfl , which is just the one caused by 
the x 2 y 2 term. The physical consequences of modality were studied in (l^] 
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that gives the dimension of the field z in the model is —2/3. However, at a = — 1 
something singular happens: The conformal dimension of z decreases to 1/16, the same 
as that of x or y, and the symmetry augments to 5*4, representing the 4-state Potts 
model Q . The potential has to include now the z 4 term, becoming 

V{x, y) = x 4 + y 4 + z 4 + xyz + w(x 2 + y 2 + z 2 ), (48) 



which was considered years ago for the 4-state Potts model in any dimension ||16|| . For 
smaller values of a, it is the term x 4 + y that has to be dropped; hence, the symmetry 
reduces to _D 4 again. 

The reader may be wondering about the bearing of the previous arguments on 
the models with p > 5. Here the field z is certainly not to be disregarded but the 
symmetry considerations still stand. In any case, the argument used for the W(%) 
potentials applies: If the term with J12 in ([41]) is omitted, the maximum unfolding 
consists of minima on spherical surfaces and the perturbation that produces it is easily 
found. The effect of J12 is enforcing the symmetry to give the pyramidal shape. A 
detailed proof would involve a construction with algebraic surfaces on the selected A 3 
diagram of dominant weights. It is straightforward and will not be expounded here. 
Finally, it is clear as well how to generalize these methods to higher values of n. 



5 Conclusions 

There are two main utilities of the lagrangian approach to 2dCFT: First, it endows 
the field algebra with further structure that constitutes a much simpler picture of 
it. Second, it allows a direct study of the phase diagram that does not need to deal 
with the complications inherent to other methods of treating perturbed 2dCFT (Bethe 
ansatz, etc). It also provides with simple objects, such as the state diagram, that can 
be directly related to statistical models known to have critical behavior described by 
those 2dCFT. It is therefore a sort of link between statistical models and 2dCFT. In 
this context, it may be of some help to people concerned with subjects that manifest 
themselves in both areas, like quantum groups. 

In this paper we have developed the lagrangian approach for the minimal models 
of W-algebras, relying on the free field construction, which directly relates them to Lie 
algebras (to be precise, Kac-Moody algebras). We have seen how all the relevant fields 
of W? n \ can be identified with monomials of the elementary fields that fit on the A n _i 
dominant- weight lattice of level 2(p — n). The main field equation follows from an OPE, 
in the Zamolodchikov's way. A new feature, worth emphasizing, is that one has to save 
in that OPE two fields instead of one to get the correct field equation. Furthermore, 
this seems to be the only requirement for consistency of the lagrangian description. We 
have obtained afterwards the Landau- Ginzburg lagrangians, noting that they posses 
the D n symmetry of these models. The potential contains the necessary information to 
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study the whole phase diagram. We have however limited ourselves to the phase struc- 
ture produced by symmetric perturbations. In particular, we have shown that there 
exists a perturbation that originates a state diagram identifyable with that defining the 
statistical models of Jimbo et al. ||, namely, the v4 ra _i dominant- weight diagram of 
level p — n. On the other hand, the phase transition between definite regimes of these 
models has already been shown to be described by the W-models. In consequence, our 
result closes the loop, showing a two-way equivalence. 

Many possibilities remain to be explored. For instance, we have confined ourselves 
to the W models that are diagonal modular invariants. Extending the methods in this 
paper to the non-diagonal ones seems quite straightforward. It would be interesting 
to see whether the resulting potentials fit in a classification (ADE or similar) agreeing 
with that known for the non-diagonal modular invariants of W models. 
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Figure Captions 



1. Level curve for a minima comfiguration with the shape of the A 2 weight diagram 
of level 3. 

2. Contour Plot of the potential V(x,y) = 0.49/ 1 2 + (/ 2 + /i)(/ -0.033), correspond- 
ing to W^y 

3. A 3 Lattice of dominant weights of level 3 with some field identifications. 



19 



